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This paper presents a method to construct variational integrators for time- 
dependent lagrangian systems. The resulting algorithms are symplectic, preserve 
the momentum map associated with a Lie group of symmetries and also describe 
the energy variation. 

Key words: variational integration algorithm, time-dependent lagrangian, symplectic 
integrator. 

1. Introduction 

There are several numerical integration methods that preserve some of the invari- 
ants of an autonomous mechanical system. In Q, T.D. Lee studies the possibility that 
time can be regarded as a bona fide dynamical variable giving a discrete time formulation 
of mechanics (see also |l^). From other point of view (integrability aspects) Veselov 
p9[ | uses a discretization of the equations of classical mechanics. Both approaches can 
be characterized as the creation of integrators based on a discretization of the varia- 
tional principle determined by a lagrangian function. These integration methods have 
usually better long term simulation properties and computational efficiency than the 
conventional ones. 

The main geometrical invariants that these integrators preserve are symplecticity, 
energy or/and momentum. Ge and Marsden have proved that a constant time step- 
ping integrator cannot preserve the symplectic form, energy and momentum, simulta- 
neously, unless it coincides with the exact solution of the initial system up to a time 
reparametrization. However, Kane, Marsden, Ortiz and West |^ show that using an 
appropriate definition of symplecticity and an adaptative time stepping it is possible to 
construct a variational integrator which is simultaneously symplectic, momentum and 
energy preserving. 
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The purpose of this paper is to extend the resuhs previously obtained for conservative 
mechanical systems to the case of time-dependent lagrangian systems following Lee's 
approach. The method is also based on a discrete variational principle. In a further 
work, we will test our algorithm in relevant examples and analyze the deep relationship 
between the continuous and discrete cases. 

From the point of view of applications to problems in mechanics (time-dependent har- 
monic oscillator, Meshchersky's equations in rocketry, for instance), and in control theory, 
the time-dependent case definitely deserves special attention. Moreover, time-dependent 
lagrangian systems appear as an indispensable tool in many economic problems JTs] , pT| . 
In fact, a typical optimization problem in modern economics deals with extremizing the 
functional D(t)U[f{t,k,k)] dt subject or not to constraints. Here, D{t) is a discount 
rate factor, U an utility function, / a consumption function and k the capital-labor ratio. 
Moreover, the time-dependent case is also important because its relation with Classical 
Field Theory, and in a future work, we will investigate the extension of these ideas 
in the framework of multisymplcctic geometry (see the preliminary work by Marsden, 
Patrick and ShkoUer 0, and ||, |ll| for variational systems subjected to constraints and 
nonholonomic mechanical systems). 

2. Time-dependent Lagrangian systems 

Time-dependent lagrangian theory is sometime studied using an homogeneous for- 
malism. If i : R X TQ — > R is a time-dependent lagrangian then the corresponding 
homogeneous lagrangian L : ro(R x Q) ^ M is defined by: 

L{t,q,i,'i) = iL{t,q,q/i) 

(we delete the submanifold t — from the whole tangent space T(M x Q).) A solution 
{t{s),q{s)) of the Euler-Lagrange equation for L verifies: 

d f dL , , . , . qis) /\ .dL, , , , , q(s) . 

d;W'^'^''^^'^'W^ 

d f^, , s , s q(s) . q(s) dL , , , , , g(s),\ , , , , , , q(s) ^ 

- [^ms), qis), |L2) - ||^(t(.), If) ) - tis)-m,qis), If) = 

Along the submanifold t = 1 these equations represent the Euler-Lagrange equations 
for L and the non-preservation of the energy function El = qdL/dq — L, respectively 
(see 1^, for instance). 

3. The time-dependent variational integrator 

A discrete time dependent lagrangian is a map M. x Q x M. x Q ^ M.. Define the 
action sum iS : (R x (3)^+^ R corresponding to the lagrangian by 

S — ^ {tk+i — tk)Ld{tk,qk,tk+i,qk+i), 

k=0 
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where qk € Q and tk gM. ior < k < N. 

We know that for any product manifold Mi x . . . x Ma, T*^^ ^ )(-^i x . . . x Ma) ~ 
T*^Mi ® ... ®T*^Ma, for all a; j G Mj, 1 < i < a. Therefore" any covector a € 
^(tci X X • • • X Ma) admits a decomposition a = ai + . . . + Ua where ai e T*.Mi, 
1 < 'i'< a. Thus, 

dLditQ,qQ,ti,qi) = DiLd{tQ,qQ,ti,qi) + D2Ld{to,qo,ti,qi) 
+D3Ld{to, qo,ti,qi)+ D4,Ld{to, qo,ti,qi) 

since T*^ (R xQxRxQ) ^ T*R ® T^^Q © T*M © T*^Q and, in addition, we 

have 

dLd{tn,qo,ti,qi) = DiLd{tn,qo,ti,qi) + D2Ld{to,qo,ti,qi) 

™ Tito,,o,tu<ii)(^ X X M X Q) = T(;„,,„)(M X 0) ©T*^^^^)(M x Q) . From now on, we 
will use both decompositions {D and D) of 1-forms on (M x Q)^. 

The discrete variational principle assures that the solutions of the discrete system 
must extremize the action sum given fixed points {to,qo) and {tN,qN)- Extremizing <S 
over tk, qk, I < k < N — 1 leads to the following system of equations: 

{tk+i — tk)D2Ld{tk,qk, tk+i,qk+i) + {tk — tk-i)D4,L(i{tk-i, qk-i,tk, qk) = (1) 

{tk+i — tk)DiLd{tk,qk, tk+i,qk+i) — Ld{tk-qk, tk+i,qk+i)dtk + 

+ {tk — tk-i)D^Ld,{tk-i,qk-i,tk,qk) + Ld{tk-i,qk-i,tk, qk)dtk = , 
or, equivalent ly, 

Di [{tk+i — tk)Ld{tk, qk, tk+i,qk+i)] + D2 [{tk — tk-i)Ld{tk-i,qk-i,tk, qk)] = . (3) 

Take Ld{tk, qk, hk, qk+i) = Ld{tk, qk, tk + hk, qk+i), where hk = tk+i - tk is the time 
step, then 

dLd _ did _ did 
dtk dtk dhk 

dLd did 



dtk+i dhk 

Therefore the second set of equations is equivalent to: 

hk-^{tk,qk, hk,qk+i) - hk-^{tk,qk,hk,qk+\) - Ld{tk,qk, hk,qk+i) + 

dtk ahk , , 

dLd - 
+/ifc-i7T7 {tk-1, qk-1, hk-1, qk) + Ld{tk-i,qk-i, hk-i,qk) = . 
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Next, we define the discrete energy function as (see Kane et al g): 
d - 

Ed{t, x,h,y) = -— [hLd{t, x, h, y)] = -Ld{t, x, h, y) - hDzLd{t, x, h, y) . 
One can also motivate this definition using the discrete variational principle (see subsec- 



tion 3.2.) 



We deduce that equation (^ is equivalent to: 

1 dljd 

T- [Ed{tk,qk,hk,qk+i) - Editk-i,qk-i,hk-i,qk)] = — 7T—{tk,qk,hk,qk+i) 

Ilk Otk 

which is a "discretization" of the well-known equation dE^/dt = —dL/dt for continuous 
time-dependent lagrangian systems, that measures the non-conservation of the energy. 
Consider the mapping 

$ : (M X Q)2 — >{Rx Qf 

defined by 

^{tk-i,qk-i,tk,qk) = {tk,qk,tk+i,qk+i), I < k < N -1 , 

where qk+i, tk+i are implicitly defined by Equations ([^). <i> will be called the discrete 
flow of the discrete lagrangian Ld- If the matrix {Di2Ld) is regular, then the discrete 
map $ exists and is uniquely defined . 

Example 1. A single particle in one-dimensional space 

Take the continuous lagrangian L(i, x, v) = "2"^^ ~ '^■^ ^^'^ ^ typical discretization 
r / , , N m f qi-qoY f qi + qo h + to 

Ld{qo,to,qi,tij ^ Y 1 _ ) 

The implicit discrete algorithm is given by: 
'qk+i-qk qk-qk-i 

m 



tk+i — tk tk — tk-i 

\dV ( qk+1 + qk tk+i+tk\ ^ dVfqk + qk-l tk + tk-i 
= "2 [{tk+i tk)^ , ^ j+{tk- tk-i)^ , 

Ak+i+tk qk+i + qk qk+i — qks ^ Ak + tk-i qk + qk-i qk — qk-i-. 

7. ' 7. ' -7 -r) ~ 7, 1 7, 1 -, 7 ) 

Z Z tk+i — tk Z Z tk — tk-i 

1 r,, dV (qk+i+qk tk+i+tk\ , , dV fqk+qk-i tk+tk-i 
{tk+i - tk)^ , j + {tk - tk-,)^ , 

where EL{t,x,v) — "^v"^ + V{t,x) is the energy function. 
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3.1. Convergence and error analysis for the discrete algorithm 

For simplicity, we suppose that Q is a vector space and choose the following dis- 
cretization of the lagrangian L -.M. x TQ M: 

T /, , \ ^Ak + tk+i Qk + Qk+i Qk+i-qk-. 

^d[tk,qk: tk+i,qk+i) — l^y Ti , , — j 

I I tk+i — tk 

In terms of the discrete velocity Vk = {qk+i — qk)/{tk+i ~ tk) and the time step hk, 

Lditk,qk,hk,Vk) = L{tk + -^,qk +Vk^,Vk) 
Then, equations (m) and (0) become 



1 

2h~k 
1 

hk 



hk^{Vk) + hk-,^{Vk-i) 
dqk dqk-l 



1 

hk 



dLd, dLd 
dvk dvk-i 



Vk^iVk) - LdiVk) - Vk-i4^iVk-i) ~ Ld{Vk-i) 

OVk UVk-l 



1 

2h^ 



hk^{Vk) + hk^i^{Vk-i) 

Otk Cltk-l 



where Vi — {ti,qi,hi,Vi). It is clear that if hk-i — > and /i^ — > then the previous 
equations converge respectively to: 



dL d fdL 
dq dt \ dq 



— and 



dt ^ dt 



As an alternative way, we consider the discretization of the homogeneous lagrangian: 

r u -t \ tk+i — tk ,tk + tk+i qk + qk+i qk+i — qks 

i~-d{tkTqk:tk+i,qk+i) — r l^y ^ , ^ , — j 



h ' 2 
where h G M+ is the time step. We then extremize 



tk+i — tk 



N-l 



N-1 



S — ^ Cd{tk,qk, tk+i,qk+i) — T ^ (ife+i — tk)L{ 



k=0 



k=0 



tk + tk+1 qk + qk+1 - qk 

2 ' 2 ' tk+1 — tk 



to derive equations ([l]) and (||) for the function 

T f, , \ j,tk + tk+i qk + qk+i qk+i-qk-. 

^d{tk,qk, tk+i,qk+i) — -t^i 7: , 7: , —) 

Z Z tk+1 — tk 

The above formalism is useful since now we can use the results of Chapter 6 of |^] to 
analyze the local truncation and global error for (|l|) and (||) as a particular case. 
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3.2. Symplecticity of the algorithm 

Now, we will show that the algorithm (|l|) is symplectic in a natural way. By similar 
arguments to those used in |^ , it is easy to prove that for solutions of the algorithm (|^) 
without imposing endpoints conditions: 

d^\solutions — {tl ~ ta)D2Ld{to,qQ, tl, 9l) + {^N — tN-l)DiLd(tN-l,qN-l,tN, In) 

+ {ti - to)DiLdito,qo,ti,qi) - Ld{to,qo,ti,qi)dto 

+ {tN — tN-l)D3Ld{tN-l,qN-l,tN, qN) + Ld{tN-l,qN-l,tN, qN)dtN 

= {ti — to)D2Ld{to, qo, ti,qi) + {t^ — tN~i)D4Ld{tN~i,qN-i,tN,qN) 
+hoDiLd{to, qo,ho, qi) + Ed{to, qo, ho, qi)dto 
—Ed{tN-i,qN-i,hj^_i,qff)dt^ . 



where hk = <fe+i — tfe. If we put 

dLd 

Energyj^ = ho-^{to,qo,ho,qi) + Ed{tQ,qQ,ho,qi) 
Energy^^ = Ed{tN-i, qN~i, hN^i^v^) 



then 

^*^| solutions 



— hoD2Ld{to,qo, ti,qi) + hN-iD4Ld{tN-i, qN-i, hN-i, qw) 
+Energyt^dto - Energyt^dtAr 



Observe the similarity of this expression with the classical transversality conditions for 
variable initial and terminal points when we are studying extremals of functionals deter- 
mined by time-dependent lagrangian systems. 

Now, fix TV = 2 and consider the 1-forms 8^^ and following a similar definition 
than in Marsden et al fl^ : 

0Zd(*o,9o,ii,9i) = (ti -to)D2Ld{to,qo,ti,qi) 

+ {ti -to)DiLdito,qo,ti,qi) - Lrf(to, go, 9i) ^^o 
0Lrf(*o,go,*i,gi) = {ti - to)D4Ld(to,qQ,ti,qi) 

+ {ti -to)D3Ldito,qo,ti,qi) + Ld{to,qo,ti,qi) dti . 



In view of the equations defining the algorithm, we deduce that 

Since ddS — 0, we have dQj^^ = -^*dQj^^, but observe that 9^^ + ^Ld = ~^o)i(i), 
that is, dQ^ + dQ^ — and, consequently, ~ Qd where fi^ = dQ^ = "'^©l • 



[Author and title] 



7 



Notice that fid ~ — dOj^, where QJ^^ is equal to 

hoDiLd{to,qo,ti,qi) - Editn,qo,ho,qi)dti 

with to + ho — ti. 

The 2-form fid is symplectic if and only if the matrix {Di2Ld) is regular. Observe 
also the similarity of the constructions of and fl^ with the Poincare-Cartan 1-form 
and 2-form, respectively, for time-dependent lagrangian systems (see p^ , p^ ) 

3.3. Momentum conservation 

The conservation of momentum for an algorithm invariant under a symmetry group 
follows in the same way as usual [po[ considering here the invariance of the function 
{tk+l — tk)Ld{tk, qk,tk+i,qk+i)- 

In particular, we suppose that {tk+i — tk)Ld{tk, qk,tk+i,qk+i) is G-invariant under 
the diagonal action of a Lie group G on M x Q. Denote by <i>g : M x Q — > R x Q the 
group action, g G G. For any element ^ in the Lie algebra g of G, we denote by ^rxq the 
fundamental vector field which generates the flow <i>cxp(s^) : M x Q ^ M x Q. Moreover 
Ckxq(^, (?) = ^M{t,q) +^Q{t,q) where ^R{t,q) G TtR and ^g(i, g) G TgQ. 

Therefore, applying the results for discrete algorithms for conservative mechanical 
systems (see for instance |20|) we obtain that the momentum map Jd : (M x Q)^ g* 
defined by 

{Jd{tk, 9fe, ife+1, ?fc+i), — {^2 [{tk+l - tk)Ld{tk,qk, ifc+i, , CRxQ(ifc+i, 9fc+i)) 

is preserved by the discrete fiow $ : (M x Q)^ — > (M x Q)^. 



4. Symplectic-Energy-Momentum preserving variational integrators 

In Kane, Marsden and Ortiz it is constructed a variational integrator for con- 
servative mechanical systems with time-step adaptation. Their algorithm is symplectic, 
energy and momentum conserving. We will show that the Kane, Marsden and Ortiz's 
algorithm is a particular case of ours when we impose that the discrete lagrangian only 
depends on the time step (hk — tk+i — tk). 

Assume that Ld{to,qo,ti,qi) = Ld{qo,qi,ti - to) = Ld{qo,qi,ho) and, therefore, 
DiLd = D2Ld, D2Ld = DiLd, and 

dLd dLd dLd 

-^('70,91,^1 - to) = ~'-^{to,qo,ti,qi) = -^(^o, go, ^i, gi) 

Then, equations (|l]) and (H) can be rewritten as 

hkDiLd{qk,qk+i,hk) + hk-iD2Ld{qk~i,qk,hk^i) = (5) 
Ed{qk,qk+i,hk) ~ Ed{qk-i,qk,hk-i) = (6) 

which are exactly the equations proposed in 

Observe that in this case the discrete lagrangian Ld is invariant by time translations, 
and the associated conservation law is precisely the energy. 
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5. Concluding remarks 

This natural variational interpretation of the results of j^] was independently obtained 
by Chen, Gou, Wu (see [||). They use Lee's approach and also derive in a variational 
way the symplectic-energy-momentum integrator for conservative lagrangian systems. 

The authors wish to thank an anonymous referee for valuable and useful comments 
and for pointing out to us reference |Q . 
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